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© Method of and device for determining spectrum parameters of a spectrum related to spectroscopic 
signals. 

© A method is proposed for the determination of 
spectrum parameters of a spectrum by quantification 
and model fitting of sampling values of a signal in 
the time domain on the basis of a model function 
comprising exponentially damped sinusoids. The 
computational efficiency of the method is very high 
with respect to known methods which are also ac- 
curate. It is proposed to start with a coarse estima- 
tion of a number of spectrum parameters, followed 
^Jby linearisation of the model function and iterative 
^estimation of the coarsely estimated parameters and 
CO further parameters, utilising a least-squares optimisa- 
tion procedure where inner product terms are 
brought into analytic form and are preferably recur- 
sively determined. 
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Method of and device for determining spectrum parameters of a spectrum related to spectroscopic 

signals. 



The invention relates to a method of determin- 
ing spectrum parameters of a spectrum related to 
spectroscopic signals, in which method the spec- 
troscopic signals derived from a substance are 
sampled in order to obtain sampling values and are 
approximated by a complex model function which 
contains the spectrum parameters and exponen- 
tially damped sinusoids, in which method there is 
also made an initial estimate, based on the sam- 
pling values, for at least one of the spectrum pa- 
rameters, the at least one and further spectrum 
parameters being accurately estimated from the 
model function by iteration with a least-squares 
optimisation procedure, and prior knowledge is in- 
troduced into the model function. 

The invention also relates to a device for deter- 
mining spectrum parameters of a spectrum related 
to spectroscopic signals, which device comprises 
means for generating the spectroscopic signals in a 
substance, sampling means for obtaining sampling 
values from the spectroscopic signals, and display 
means for displaying the spectrum, and also com- 
prises programmed means for making an initial 
estimate for at least one spectrum parameter, 
which programmed means also comprise a model 
function of exponentially damped sinusoids and are 
suitable for storing prior knowledge of the spec- 
troscopic signals and for executing a least-squares 
optimisation procedure on the basis of the sam- 
pling values, the programmed means furthermore 
comprising iteration means for accurately estimat- 
ing, the at least one and further spectrum param- 
eters by means of the least-squares optimisation 
procedure. 

A method of this kind is inter aiia suitable for 
signals which are mainly exponentially damped, for 
example for signals obtained during a magnetic 
resonance experiment from an entire body as the 
substance as well as from a part of the body. The 
method can also be used, for examle for X-ray 
spectroscopy or FT infrared spectroscopy. 

A method of this kind is disclosed in an article 
by J.W.C. van der Veen et al. p "Accurate Quan- 
tification of in Vivo 31 p NMR Signals Using the 
Variable Projection Method and Prior Knowledge", 
Magnetic Resonance in Medicine, Vol. 6, No. 1, 
January 1988, pp.92-98. The spectrum parameters 
of a spectrum related to spectroscopic signals 
(signals obtained during an NMR experiment) are 
derived via the sampling values (represented in a 
time domain, which is In contrast with the spectrum 
which is represented in a frequency domain), di- 
rectly in the time domain. A least-squares op- 
timisation procedure is used for fitting the sampling 



values to a model function of spectroscopic sig- 
nals, in which procedure prior knowledge concern- 
ing spectral components can be introduced. The 
method can in principle be used for an arbitrary 

5 model function. In said article, for example a model 
function is assumed which contains exponentially 
damped sinusoids. By operation in the time do- 
main, arbitrary parts of the spectroscopic signals 
can be omitted without giving rise to serious fitting 

70 problems. For example, when in the case of NMR 
the substance contains immobile nuclei and mobile 
nuclei, the sampling values (of a quickly exponen- 
tially decaying spectroscopic signal) originating 
from the immobile nuclei can be simply omitted. 

15 The sampling values of the tail of an exponentially 
decaying spectroscopic signal can also be omitted. 
Thus, convolution effects associated with a trans- 
formation to the frequency domain in conjunction 
with weighting, such as base line (wide background 

20 in the spectrum due to quickly decaying signals 
from the immobile nuclei) and line-shaped distor- 
tions (due to phase shifts) can be avoided. Fitting 
is performed partly non-iteratively in order to obtain 
an initial estimate for at least one of the spectrum 

25 parameters used as initial values for further itera- 
tive fitting. The non-iterative fitting is comparatively 
fast and the further iterative fitting enables the 
method to be used successfully also in the case of 
a poor signal-to-noise ratio. For parameters which 

30 occur linearly in the model function, such as com- 
plex amplitudes, no initial values are required, re- 
sulting, in a faster method. Said article describes 
the method inter alia for a model function of ex- 
ponentially damped sinusoids in which parameters 

35 such as amplitude, damping factor and frequency 
occur. Using the least-squares optimisation proce- 
dure, the sampling values are iteratively fitted as 
well as possible to the model function, initial values 
being required for the damping factor and fre- 

40 quencies which occur as non-linear parameters in 
the model function. The prior knowledge may com- 
prise a time shift of an instant of a first sample with 
respect to a time origin and predetermined phases. 
The least-squares optimisation procedure can be 

45 written in the form of a matrix. This is represented 
by the formule 4 on page 94 of said article, show- 
ing a least-squares solution for the amplitude. All 
terms in the inverse matrix occurring therein are 
successively calculated, like the terms in the prod- 
so uct of the Hermitic matrix with the data vector. In 
the matrix form a large number of product terms 
occur, so that a substantial amount of computation 
time is required for executing the method by 
means of a computer; this aspect becomes more 
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significant in the case of a comparatively large 
number of sampling values. 

It is an object of the invention to provide a 
method which, when the model function contains 
exponentially damped sinusoids, requires substan- 
tially less computation time when executed by 
means of a computer. 

To achieve this, a method in accordance with 
the invention is characterized in that prior to the 
iteration the model function is linearised with re- 
spect to the at least one spectrum parameter, after 
which inner product terms occurring in the least- 
squares procedure are internally converted into an 
analytic form which is used during iteration. The 
number of columns of the matrix is increased with 
respect to the matrix in said article due to residual 
parameters which remain after linearisation, for ex- 
ample frequency and damping. The invention is 
based on the idea to convert first the inner product 
terms into analytic form in order to solve the least- 
squares problem in the case of exponentially 
damped sinusoids as the model function, rather 
than calculate, like in the known method, all matrix 
products one term after the other for all possible 
summations by means of row and column mul- 
tiplications. The computation time is thus substan- 
tially reduced. In comparison with the known meth- 
od, the method in accordance with the invention is 
30 times faster in the case of 17 sinusoids and 512 
sampling values. Random tests have shown that 
the Cramer-Rao lower limit for estimating the pa- 
rameters is at least substantially reached, so that 
the method produces reliable parameters. The 
computation time increases in proportion to the first 
power of the number of sampling values, assuming 
that the other parameters, for example the number 
of sinusoids, remain the same. 

A version of a method in accordance with the 
invention is characterized in that terms of the ana- 
lytic form are recursively determined. Thus, the 
method can be simply implemented in a computer. 

A version of a method in accordance with the 
invention is characterized in that the model function 
contains a polynomial whereby the exponentially 
damped sinusoids are multiplied. A sum of ex- 
ponentially damped sinusoids does not always con- 
stitute a good model function in all circumstances. 
By multiplication by a polynomial, a model function 
is obtained which is adapted to circumstances. The 
inner product terms retain their comparatively sim- 
ple form. 

A version of a method in accordance with the 
invention is characterized in that the initial estimate 
of the spectrum parameters is made by means of a 
non-iterative estimation procedure. Initial values 
can. thus be determined for parameters which oc- 
cur non-linearly in the model function. For a non- 
iterative estimation procedure HSVD (Hankel matrix 



singular value decomposition) can be used as de- 
scribed in "Improved Algorithm for Noniterative 
Time-Domain Model Fitting to Exponentially 
Damped Magnetic Resonance Signals", H. Bark- 

5 huijsen et al, Journal of Magnetic Resonance 73, 
pp. 553-557, 1987. 

A preferred version of a method in accordance 
with the invention is characterized in that the initial 
estimate of the spectrum parameters is made on 

w the basis of a coarse spectrum obtained from the 
sampling values by Fourier transformation. In the 
case of a comparatively poor signal-to-noise ratio, 
first a filtering operation can be performed. Initial 
values of the parameters which occur non-linearly 

75 in the model function can be determined from the 
coarse spectrum. These parameters, frequency 
and damping, can be determined directly from a 
display screen on which the coarse spectrum is 
displayed by picking up, using a so-called mouse, 

20 at peak locations in the coarse spectrum the peak 
locations themselves (frequencies) and the peak 
width at half peak height (a measure of the damp- 
ing) for supply to the programmed means in a 
computer. When the coarse spectrum contains ex- 

25 cessive phase distortion, first a phase correction 
method incorporated in the programmed means 
can be executed, possibly after filtering (in the 
case of a comparatively poor signal-to-noise ratio), 
for example as described in an article by C.H. 

30 Sotak et al in Journal of Magnetic Resonance, Vol. 
57, pp. 453-462, 1984. The phase correction of the 
coarse spectrum and the determination of the initial 
values can also be performed by means of a 
method as described in the noh-prepublished Neth- 

35 erlands Patent Application No. 8702701 (PHN 
12325), "Werkwijze en inheriting voor automatische 
facecorrectie van complexe NMR spectra", by J. 
van Vaals. 

The invention will be described in detail 
40 hereinafter with reference to a drawing, therein 

Fig. 1 diagramrnatically shows a device in 
accordance with the invention, 

Fig. 2 shows, in the form of a matrix, a set of 
linear equations in accordance with the invention, 
45 to be used for each iteration step during the iter- 
ation process, 

Fig. 3 shows the analytic form to be substi- 
tuted, in accordance with the invention, in the set of 
linear equations given in Fig. 2, 
so Fig. 4 shows a spectrum obtained by means 

of the method in accordance with the invention, 
and 

Fig. 5 shows a table containing parameters 
of the spectrum of Fig. 4. 

55 

Fig. 1 diagramrnatically shows a device 1 in 
accordance with the invention. Specifically, a mag- 
netic resonance device for obtaining spectroscopic 
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signals is shown by way of example. The device 1 
comprises magnet coils 3 for generating a steady, 
uniform magnetic field B 0 which are arranged with- 
in a shielded space 2. Within the magnet coils 3 
there may be arranged a substance. The substance 
is, for example a body containing solid and liquid 
components. During operation of the device, the 
substance being arranged within the magnet coils 
3. a slight excess of nuclear spins (of nuclei having 
a magnetic moment) will be oriented in the same 
direction as the field B 0 . From a macroscopic point 
of view this is to be considered as a magnetisation 
M, i.e. anequilibrium magnetisation. The magnet 
coils 3 may be resistance magnets or supercon- 
ducting magnets; the device then comprises a d.c. 
power supply source 4 for powering the magnet 
coils 3. When the magnet coils 3 are constructed 
as permanent magnets, the d.c. power supply 
source 4 will be omitted. The device 1 furthermore 
comprises means 5 for generating spectroscopic 
signals, which means include a reference generator 
6 for generating a reference signal for supply to an 
r.f. modulator 7. When activated, the r.f. modulator 
7, supplied an r.f. electromagnetic pulse to an r.f. 
power amplifier 8 which amplifies the r.f. elec- 
tromagnetic pulse for supply to a transmitter coil 9 
which is arranged near the magnet coils 3. When 
the substance arranged within the magnet coils is 
irradiated, the equilibrium of the magnetisation M is 
disturbed and after irradiation there will be gen- 
erated a magnetic resonance signal when the mag- 
netisation returns to the state of equilibrium, a so- 
called FID signal (Free Induction Decay). The reso- 
nance signal contains spectroscopic information: 
the signal will be referred to hereinafter as a spec- 
troscopic signal. The device 1 also comprises 
means 10 for receiving the spectroscopic signals, 
including a receiver coil 11 for receiving the spec- 
troscopic signals which is connected to a detector 

12 for the detection of the modulated spectroscopic 
signal, for example by quadrature detection. The 
detector 12, being coupled to the reference gener- 
ator 6, comprises low-pass filters and analog-to- 
digital converters as the sampling means for sam- 
pling the spectroscopic signals received and de- 
tected. The device 1 also comprises control means 

13 which serve inter alia for control and timing of 
the r.f. generator 7. The device furthermore com- 
prises processing means 14 for processing digitis- 
ed spectroscopic signals 15 (sampling values). The 
processing means 14 are coupled to the control 
means 13. The processing means 14 comprise 
programmed means 16 inter alia for making an 
initial estimate for at least one spectrum parameter 
of a spectrum related to the spectroscopic signals 
in order to perform a least-squares optimisation 
procedure on the basis of the sampling values. The 
programmed means 16 furthermore contain a 



model function of exponentially damped sinusoids, 
a non-iterative estimation procedure, iteration 
means and Fourier transformation means (not 
shown in detail). The processing means 14 also 

5 comprise a RAM memory 17 which is coupled to 
the programmed means 16. The device 1 further- 
more comprises display means 18, inter alia for the 
display of a spectrum of the spectroscopic signals. 
A cursor control device 19 is coupled to the display 

w means 18 and the processing means 14. The cur- 
sor control device may be a so-called mouse or 
tracker ball; as is known, the mouse is formed by a 
carriage provided with a motion sensor which out- 
put signals in response to motion of the mouse. 

rs The signals of the motion sensor are applied to the 
processing means 14 which move a cursor as an 
indicator across a display screen of the display 
device in dependence of the signals of the motion 
sensor. When a given screen position is reached, 

20 the screen position can be applied to the process- 
ing means by operation of a control member on the 
mouse. Inter alia for spectroscopy of a part of the 
substance, so-called volume-selective spectro- 
scopy, the device 1 comprises gradient coils 20 

25 which are controlled by a gradient coil control 
device 21, The gradient coil control device 21 is 
controlled by the control means 13. The gradient 
coils 20 serve to generate field gradients which are 
superposed on the field Bo. Generally, there can be 

30 generated three gradients Gx, Gy. and Gz whose 
field direction coincides with the' direction of the 
steady, uniform magnetic field Bo and the respec- 
tive gradient directions of which extend mutually 
perpendicularly. For a more detailed description of 

35 a magnetic resonance device, reference is made to 
the handbook "Practical NMR Imaging", by M.A. 
Foster and J.M.S. Hutchinson, IRL Press, 1987, 
ISBN 1-85221-011-7, Chapter 1, pages 1-48 which 
describes the general principles of NMR and 

40 shows an NMR device. In the case of spectroscopy 
of the substance or part of the substance, the aim 
is to obtain information from the substance in the 
form of a spectrum (in the frequency domain). For 
NMR, where nuclear spins are excited in the 

45 steady field Bo, the relation <u =gamma.(1-a).Bo 
holds for excitation of nuclear spins, where w is a 
resonant frequency of the nuclear spins, gamma is 
a gyromagnetic ratio, and a is a so-called shielding 
constant of a nucleus. Spin resonance will occur 

so when the resonant frequency is applied. From a 
macroscopic point of view this implies, vis-a-vis a 
stationary observer, a precessional motion of the 
magnetisation M around the field Bo at the reso- 
nant frequency. When the substance contains 

55 nuclei having a different o, the resonant frequency 
will also be different for these nuclei. This effect is 
referred to as chemical shift and becomes manifest 
as different resonant peaks in a spectrum obtained 
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from the resonance signal. In the time domain the 
resonance signal can in principle be described as a 
superposition of exponentially damped sinusoids. 
The chemical shifts are due to the fact that nuclei 
effectively experience a different magnetic field 
due to a different shielding under the influence of 
the presence of neighbouring electrons. The neigh- 
bouring eiectrons shield the nuclei to a given de- 
gree from the field Bo. Thus, chemical components 
can be identified in a substance by means of NMR 
spectroscopy. Depending on a selected frequency 
of the r.f. magnetic pulse, a given spectrum will be 
obtained. For example, a phosphor spectrum ( 31 P), 
a carbon spectrum ( 13 C) or a proton spectrum ('H) 
can be obtained. The surface area below the reso- 
nance peaks in the spectrum is directly propor- 
tional to the number of nuclei making a contribution 
for a given frequency. The surface area corre- 
sponds to the amplitude in the time domain. In 
addition to chemical shifts, another interaction also 
has an effect on the spectrum, the so-called spin- 
spin coupling which, which causes the splitting off 
of a resonance line which is usually small in com- 
parison with the chemical shift. For a more detailed 
description of NMR spectroscopy reference is 
made to pages 249-274 of said handbook by Fos- 
ter and Hutchinson. 

Fig. 2 shows a set of linear equations (2A) to 
be used in accordance with the invention, for each 
iteration step during the iteration process. After an 
initial estimate of at least one spectrum parameter, 
frequency and damping of peaks occurring in the 
spectrum in the present example, the parameters 
frequency and damping and also amplitude are 
accurately estimated in a number of iteration steps 
on the basis of the set (2A). Therein, the column 
vector chelates to the amplitudes of the peaks of 
the spectrum and Re c~~ represents the residual 
parameters Ap of the damping p and Im ^repre- 
sents the residual parameters Anu of the frequency 
nu, M1 to M6 being matrices in the partitioned form 
of the extended Fisher data matrix M of the least- 
squares problem, extended with columns for the 
residual parameters Ap and Anu for ail peaks, and 
v~t to v~a a representation of the sampling values. 
The set (2A) is formed by linearisation of the least- 
squares problem. The procedure is based on a 
model function for the magnetic resonance signal 
x n of the form: 

x n = E^expKpk + 27ri.nu k )(n + 5) + i<f> k ] (1 ) 
where x n is a model function of the signal x n , Ck, 
nu k and p k are parameters of the spectrum van x n , 
being amplitude, frequency and damping, respec- 
tively, for k = 1 K, in which K is the number of 

peaks in the model function. Furthermore, i = (-1) 1/2 . 
First an initial estimate is made for p k and nu k . In 
the working point (po, nuo) obtained for the param- 
eters p and nu, for each k x n is linearised, produc- 



ing terms of the form: 

[c +Jn + 5)c'].exp[(po + 27ri.nu 0 )(n + 5) + ip Q ] (2) 
in x n .p= po+Ap and nu = nuo + Anu. Therein, n 
represents a sampling time of the resonance sijj- 

5 nal, n = 0, 1 N-1 for N sampling operations at. 

instants t n ,' & denotes a time shift caused by in- 
strumental limitations and/or the deliberate omis- 
sion of sampling values, fa denotes a phase shift. 5 
and fa can be introduced into the model x n as prior 

70 knowledge. Furthermore, after linearisation Re(c )- 
= cAp and Im(c') = 2 -n-cAnu, which can be jsimply 
demonstrated. Formules (2B) to (2G) follow from 
further elaboration of the least-squares problem 
after linearisation. In the formule (2A) to (2G), ~ is a 

75 transposition and ' is a complex conjugated opera- 
tion. f k (t n ) is the k th model function, The information 
matrix is dimensioned to 3Kx3K. Furthermore: 
fk(W ■ exp[( Pk + 277i.nu k )(n + 5) + i<*> k ] (3) 
After each iteration step the current values of p k 

20 and nuk are replaced, utilising the current values of 
c k '. When c~' is sufficiently small, the iteration 
stops. A spectrum is obtained by transforming the 
model obtained after iteration to the frequency do- 
main, for example by means of the Fourier trans- 

25 formation means. The spectrum can also be cal- 
culated directly from the spectrum parameters ob- 
tained, using known Lorentz formules, so that any 
artefacts due to transformation are avoided. The 
spectrum is displayed by means of the display 

30 means 18. The spectrum parameters- can also be 
represented in the form of a table; in that case the 
Fourier transformation may be omitted. Random 
tests have demonstrated that the so-called Cramer- 
Rao lower limit is at least substantially reached. For 

35 a further description of the Cramer-Rao lower limit 
in relation to the Fisher data matrix reference is 
made to "Handbook of Measurement Science", 
Vol. 1, P.H. Sydenham, 1982, Wiley, pp. 335-339. 
For a description of the generally known Fourier 

40 transformation means reference is made to the 
handbook "An Introduction to the Analysis and 
Processing of Signals" P.A. Lynn, McMillan 1983, 
pp. 57-64. 

Fig. 3 shows the analytic form for the inner 
45 product terms to be substituted in the set of linear 
equations shown in Rg. 2. By elaborating in the 
formules (2B) to (2G) the product fk(tn)f k '*<t n ), multi- 
plied by a power of t n , for f k (t n ) in accordance with 
(3), terms according to (3A) in Rg. 3 are produced 
so in as far as the part thereof which is dependent on 
the sampling instant is concerned; therein, 1 in- 
dicates the formules corresponding to the formules 
(2B), (2C) and (2D) in Rg. 2. In formule (3A): 
2 = exp[pk+pk'+27ri(nu k -nuk / )] (4) 
55 Formules (3B) and (3C) show recursive relations for 
the respective terms corresponding to l = 0, 1 and 
2. When z approximates 1, an inaccuracy may 
occur when use is made of the recursive formules 
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(3B) and (3C). By introducing the variable eps = z- 
1, this problem an be circumvented. (3D) to (3F) 
show formules for use in such a case. The com- 
putation time required will then slightly increase 
again, notably in the case of a large N, because of 
the combinations occurring therein, for example 
when N exceeds n + 1. Without substantial loss of 
accuracy truncation can take place for a suitable 
value of n. so that the computation efficiency is 
increased again. For example, in the case of 512 
sampling values it will suffice to truncate the sum- 
ming operations for, for example n = 10. 

Fig. 4 shows a spectrum obtained by means of 
the method in accordance with the invention for in 
vivo * 3 C NMR spectroscopy of the liver of a rat 
The number of complex sampling operations 
amounted to 512 with a sampling interval of 50 us. 
The frequency scale is given in kHz relative to the 
irradiation frequency (the chemical shift appears as 
a spectrum). In Fig, 4 the reference 4A is the 
cosine FFT (Fast Fourier Transform) of the sam- 
pling values { a window is used in order to avoid 
artefacts), the reference 4B is a spectrum obtained 
by means of a method other than the method in 
accordance with the invention, and the reference 
4C is the cosine FFT of the model which has been 
fitted in the time domain by means of the method 
in accordance with the invention and in which the 
prior knowledge concerning the phases has been 
introduced. In accordance with the invention, one 
operation step required 15 seconds, * using the 
same device, whereas the known method required 
a period of time which is 30 times longer. The 
references p1 tp p17 denote resonance peaks. The 
spectrum shown in Fig. 4C was obtained by means 
of 8 iteration steps. 

Fig. 5 shows a table containing the parameters 
of the spectrum 4c in Rg. 4. The number of 
exponentially damped sinusoids was initially set at 
25; it was found that 17 sinusoids were of spec- 
troscopic relevance. For each resonance peak 
there are shown parameters obtained using a 
method other than the method in accordance with 
the invention, in this case by means of HSVD 
without utilising prior knowledge and therebelow 
the parameters obtained utilising the method in 
accordance with the invention with introduction of 
prior knowledge. The errors indicated amount to 
twice the standard deviation. 

Many alternatives are feasible for those skilled 
in the art without departing from the scope of the 
invention. For example, in the case of NMR many 
pulse sequences can be used for spectroscopy in 
order to obtain resonance signals, both volume 
selective as well as non-volume selective. For a 
volume-selective pulse sequence use can be made 
of the SPARS method, as disclosed in Magnetic 
Resonance Imaging, Vol. 4, pp.237-239, 1986. 



Even more prior knowledge can be used, for exam- 
ple information concerning multiplets in the spec- 
trum. For a given signal-to-noise ratio for in vivo 
NMR, the use of prior knowledge will often be 
5 required in order to achieve suitable quantification. 



Claims 

w 1. A method of determining spectrum param- 
eters of a spectrum related to spectroscopic sig- 
nals, in which method the spectroscopic signals 
derived from a substance are sampled in order to 
obtain sampling values and are approximated by a 

;s complex model function which contains the spec- 
trum parameters and exponentially damped 
sinusoids, in which method there is also made an 
initial estimate, based on the sampling values, for 
at least one of the spectrum parameters, the at 

20 least one and further spectrum parameters being 
accurately estimated from the model function by 
iteration with a least-squares optimisation proce- 
dure, and prior knowledge is introduced into the 
model function, characterized in that prior to the 

25 iteration the model function is linearised with re- 
spect to the at least on© spectrum parameter, after 
which product terms occurring in the least-squares 
procedure are inner converted into an analytic form 
which is used during iteration. 

30 2. A method as claimed in Claim 1 , character- 
ized in that terms of the analytic form are recur- 
sively determined. 

3. A method as claimed in Claim 1 or 2, 
characterized in that the model function contains a 

35 polynomial whereby the exponentially damped 
sinusoids are multiplied. 

4. A method as claimed in Claim 1, 2 or 3, 
characterized in that the initial estimate of the 
spectrum parameters is made by means of a non- 
40 iterative estimation procedure. 

5. A method as claimed in Claim 1, 2 or 3, 
. characterized in that the initial estimate of the 

spectrum parameters is made on the basis of a 
coarse spectrum obtained from the sampling val- 
45 ues by Fourier transformation. 

6. A device for determining spectrum param- 
eters of a spectrum related to spectroscopic sig- 
nals, which device comprises means for generating 
the spectroscopic signals in a substance, sampling 

50 means for obtaining sampling values from the 
spectroscopic signals, and display means for dis- 
playing the spectrum, and also comprises pro- 
grammed means for making an initial estimate for 
at least one spectrum parameter, which pro- 

55 grammed means also comprise a model function of 
exponentially damped sinusoids and are suitable 
for storing prior knowledge of the spectroscopic 
signals and for executing a least-squares optimisa- 
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tion procedure on the basis of the sampling values, 
the programmed means furthermore comprising it- 
eration means for accurately estimating, using iter- 
ation, the at least one and further spectrum param- 
eters by means of the least-squares optimisation s 
procedure, characterized in that the iteration means 
contain the model function linearised with respect 
to the at least one spectrum parameter, product 
terms being inner converted into an analytic form 
in the least-squares optimisation procedure, which w 
analytic form serves for use in the iteration means. 

7. A device as claimed in Claim 6, character- 
ized in that the programmed means also contain a 
non-iterative estimation procedure for making an 
initial estimate of the at least one spectrum param- is 
eter. 

8. A device as claimed in Claim 6 which com- 
prises Fourier transformation means, characterized 
in that the programmed means are also suitable for 
determining, in conjunction with the Fourier trans- 20 
formation means, a coarse spectrum from the sam- 
pling values and for displaying the coarse spec- 
trum by way of the display means, the device 
comprising a cursor control device which is coup- 
led to the display device in order to pick up param- 25 
eters from the coarse spectrum as an initial es- 
timate for the spectrum parameters. 
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